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PURPOSE 


The  tests  performed  prior  to  this  test  run  established  the  technical 
and  economic  feasibil ity  of  efficient  fluoride  removal  via  alumina  column 
adsorption.  Th-is  followed  with  second  generation  research  effort  that  has 
or  is  being  completed  on  aluminas  column  adsorption  modes,  but-the 
d^riis|l^|)]ij|^sa  predictive  mathematical  model  for  alumina  columns  had 
not-been  made,  These  tests  were  made  to  find  a  predictive  mathematical 
model  for  alumina  columns.  This  predictive  model  uses  a  psuedo-kinetic 
rates  constant,- rate  expression  which  enables  alumina  column  behavior  to 
be  predicted  for  selected  aluminas.  The  predicted  behavior  (theoretical) 

j  ■  /  I 

can-  "Wien  be  compared  to  the  actual  behavior  to  determine  if  the  actual  - 
behavior  follows  the  predictive  behavior.  Significant  deviations  would 
indicate  either  an  alumina  column  plant  malfunction  and/or  mi  sad justed;^ 
£hus  an  early  diagnosis  (warning)  could  be  made  of  plant  failure  and 
subsequent  fix  would  greatly  shorten  downtime  and  improve  the  overall 
plarcfc  efficiency. 


DETERMINATION  OF  KINETIC  RATE  CONSTANTS  IN  ALUMINA  ADSORPTION  FOR 

SELECTED  ALUMINAS 

INTRODUCTION 

The  need  to  find  kinetic  rate  constants  for  selected  aluminas  was 
.generated  from  the  alumina  column  test-results  at  Rocky  Mountain  Arsenal. 
The  problem  that  had vnot  been  addressed  was  finding  a  means  of  calculating 
the  initial  fluoride  removal  by  aluminas  before  the  linear  portion  of  the 
fluoride  removal  versus  gallonage  curve  had  been  reached.  It  is  important 
especially  with  alumina  adsorption  columns  to  define  the  above  mentioned 
region.  The'reason  is  improper  switching  of  alumina  adsorption  columns* 
.could  result  in  a  loss  in  efficiency.  Information  about  the  region  before 
linear  fluoride  removal  was  reached  would  reduce  the  loss  inefficiency, 
accordingly,  test-runs  were  scheduled.  Since  it  was  not  known  which  of 
the  available  aluminas  would  be  used  in  the  final  design  configuration* 
six  aluminas  were  chosen  for  a  pseudo-kinetic**  rate  constant  test  runs. 

The  rate  adsorption  constants  and  rate  expressions  derived  would  allow  the 
prelinear  region  to  be  approximated  closely.  The  proper  switching  of 
alumina  columns  could  then  be  made  based  on  the  values  for  the  items 
above.  The  experimental  and  mathematical  methodology  are  given  in  the 
procedure  section. 

*  It  is  anticipated  two  or  mole  alumina  columns  will  be  needed  to  allow  a 
regeneration  cycle  to  be  performed,  especially  if  pH  adjustment  used. 

**  pseudo- kinetic  or  apparent  rate;  the  rate  of  adsorption  is  actually  much 


faster. 


EXPERIMENTAL  PROCEDURE 

The  Items  needed  to  determine  a  kinetic  rate  adsorption  constant 

and  a  rate  expression  are  the  initial  concentration  (C/^)  and  the 

concentrations  (C^)  at  set  elapsed  times  (t).  The  generalized  rate 

expression  is  +rA  =  -dC&,  that  is  +rf  =  -dCa  =  kf  (C.4)  with  rf  the 

dt  dt 

pseudo-kinetic  rate  of  adsorption;  C4  the  concentration  at  time  t. 

The  intergrated  rate  expression  is  then  ^*dCa/f(Cg)  =  -k^"dt  with  k 
the  kinetic  rate  constant  and  t  the  elapsed  time.  Temperaure,  pressure 
and  alumina  dosage  were  maintained  constant  with  time  of  contact  being 
varied.  The  alumina  dosage  was  selected  so  that  inter-ionic  competition 
would  be  minimized  between  fluoride  ions. 

To  find  C:a5  and  C^at  t  each  rate  test  was  done  in  the  following 
manner: 

(1)  Eight  jars  received  5  grams  of  the  alumina  being  tested  with 
a  ninth  jar  set  aside  for  a  control; 

(2)  The  nine  jars  were  filled  with  800  ml  of  Calgon  plant  effluent 
water,  placed  in  a  constant  temperature  bath  with  agitators  started; 

(3)  Each  hour  a  jar  was  removed  and  immediately  filtered  to  'stop' 
the  adsorption  action  of  the  alumina,  the  control  (blank)  was  filtered 
with  the  eight  hour  elapsed  time  sample; 

(4)  Fluoride  concentration  and  pH  measurements  were  found  on  the 
samples  and  recorded  along  with  elapsed  time  for  the  sample. 


MATHEMATICAL  PROCEDURE 

ttefininij  as  the  apparent  rate  of  disappearance  of  fluoride  ion, 

Cp  as  the  concentration  of  fluoride  ion  in  solution  and  t  as  the  time, 

the  apparent  kinetic  rate  is  given  by  -dCp  =  +rjp  for  the  reduction 

_  . .  dt  .  . 

of  fluoride  ion;  this  can  also  be  written  -dCp  =  kf-(Cp).  The 

dt  - 

equations  given  above  are  the  generalized  rate  expressions, however, 

in  order  to  obtain  the  integrated  rate  expression  for  predicting  the 

adsorption  of  fluoride  ions  from  solution,  f(Cp)  must  be  found. 

Rearranging  the  pseudo-kinetic  rate  expression  gives  -dCF  =  kdt  so  that 

fW) 

the  integrated  rate  expression  can  be  found,  when  f  (Cp)  known,  by 

N 

integration.  One  method^is  to  assume  f(Cp)  has  the  form  Cp  where 
N  *  O,  T/G,  1,  G  (G  =  real  number)  which  limits  the  form  of  the 
integrated  pseudo-kinetic  rate  expression  to  four  forms.  The  four 
Integrated  pseudo-kinetic  rate  expressions  are  then: 

0)  ca0  =  initial  concentration  and  Ca  =  concentration  at  time  t. 

(C3'Ca0)  =  -kt  +  Ic  (Ic  -  constant  of  integration), 

(2)  G  (Ca  -  Cao  1/G)  =  -kt  +  Ic  ,  .......... 

(31.  1n.(Ca/Cao)  =  -kt  +  Ic  and  .  ... 

,,..(4) .  ys.(ca-G+1  -  ca0‘G+1)  =  -kt  +  ic 

Now  Cao  (the  initial  concentration)  is  a  constant  and  the  Cao  terms 
^q-.Cao  are  therefore  constant  and  can  be  represented  by  Kc  (a 
constant  term) . 


- - - - - H ““  1 

■  | 

i 

The  four  pseudo-kinetic  rate  expressions  then  become  (after  substitution 
and  rearrangement) : 

(1)  Ca  =  Kc  +  Ic  -  kt, 

(2)  G  Ca1/6  •  Kc  +  Ic  -  kt, 

(3)  In  +  Iq  “  kt,  and 

(4)  1  c  -G  +  1  =  Kc  +  Ic  -  kt.  j 

G  1 

Further  combining  the  constants  Kc  and  Ic  to  give  a  new  constant  Jc  -  Kc  +  Ic 
renders  the  four  pseudo-kinetic  rate  expressions  in  the  final  forms: 

0)  Ca  +  Jc  -  kt, 

(2)  G  Ca  1/G  =  Kc  -  kf, 

(3)  In  Ca  =  Jc  -  kt,  and 

(4)  1  C  -G  +  1  =  Jc'-  kt. 

G 

The  problem  now  becomes  one  of  selecting  one  of  the  above  equations  which 
represents  the  best  fit  to  the  Ca  versus  t  plot.  To  select  the  best  fit, 
the  tabulized  data  was  taken  to  MISO  and  analyzed  via  a  modified  simple 
regression  program.  This  modified  simple  regression  program  calculated 
the  correlation  coefficient*,  residuals,  plotted  the  equation  selected 
from  the  regression  analysis,  gave  the  95%  confidence  bands  for  the  linear 
case,  and  gave  the  k  G  and  Ic  values.  Since  the  aluminas  are  constrained 
to  follow  the  same  rate  mechanism  (by  physical  chemistry  laws),  it 
was  expected  one  f  (Ca)  would  be  found  which  would  yield  a  high 
correlation  coefficient  for  all  the  aluminas. 

% 

!  * 

*  Correlation  coefficient  (R2)  =  1.00  -  direct  correlation  (X  =  y), 
excellent  fit;  R2  =  -1.00  -  inverse  correlation  (1/x  =  y),  excellent  equation 
fit;  R2  =  0,0  -  no  correlation,  R2  between  0.0  and  _+  1.00  -  partial 
correlation  with  stronger  correlation  as  R2  _ ^  +  1.00. 


r  . 


DATA  AND  RESULTS  OF  REGRESSION  ANALYSIS 

The  concentrations  of  fluoride  ion  found  for  each  elapsed  time  for  the 

six  aluminas  tested  are  given  in  Table  I  below: 

TABLE  I :  CONCENTRATION  OF  FLUORIDE  ION  VERSUS  ELAPSED  TIME  FOR  SIX  SELECTED 
ALUMINAS 


ALUMINA 
Elapsed  Time 
Hours 


MCB  Alcoa  Alcoa  Kaiser  Kaiser  Kaiser 

(crushed)  F-l  F-l  A-300  A-201  8  Mesh 

(Spheres)  (ReaeneratedXVirgin) (Ungraded) (Spheres)  (Spheres), 

;  CONCENTRATION' OF  FLUORIDE  ION  (PPM) 


1  * 

2.69 

3.08 

2.52 

1.48 

3.35 

3.12 

2 

2.40 

2.92 

1.97 

0.978 

3.26 

2.94 

3 

2.00 

2.83 

1.84 

0.780 

3.18 

2.86 

4 

1.68 

2.75 

1.49 

0.778 

3.12 

2.48 

5 

1,46 

2.73 

1.43 

0.766 

3.04 

2.24 

6 

1.02 

2.70 

1.41 

0.729 

3.02 

2.17 

7 

1,00 

2.60 

1.28 

0.608 

3.00 

2.15 

8 

0.907 

2.59 

1.19 

0.560 

2.95 

1.69 

8  (Blank) 

* 

(initial  Concen) 

3.60 

3.71 

3.62 

3.60 

3.60 

3.74 

»*v<j 


The  regression  analysis  was  initially  performed  with  Cao  (the  initial 
concentration  of  fluoride  ion)  included.  Two  conclusions  were  drawn  from 
an  analysis  of  the  regression  calculations  :  * 

(1)  overall  f  (Ca)  could  take  any  of  the  four  forms  offered  previously 

(2)  the  spheres  all  had  Ca  =  Ke  exp  (-kt)  as  the  1st  or  2nd  best 

% 

correlation. 

The  initial  regression  results  with  Cao  included  are  given  in  Figures 
1-6.  The  plots  of  Ca  vs.  Ke  exp  kt  for  the  six  aluminas  given  in  Figures 

7-12. 


i  'i 


r  .r 


When  It  comes  to  adsorption,  the  aluminas  are  chemically  the  same  but 
physically  the  aluminas  are  different.  The  aluminas  fall  into  two  physical 
groups,  spheres  and  granulated  powders.  The  great  difference  in  immediate 
surface  area  was  postulated  as  the  primary  reason  the  spheres  followed  a  In 
concentration  (or  exponential  kt)  form  while  the  powders  seemed  to  follow 
poorer  formulations.  In  order  to  determine  the  behavior  of  the  aluminas 
without  the  effect  of  large  initially  available  surface  area  versus 
small  initially  available  surface  area,  which  occurs  during  the  start-up  of 
the  tests,  the  zeroeth  and  first  hour  concentrations  were  deleted  from  the 
granulated  powder  regression  analyses.  Also  to  show  that  the  small  surface 
Of  the  spheres  had  very  little  effect  on  the  correlation  coefficient,  thus 
the  fit  of  the  equation,  the  zeroeth  hour  (initial  concentration)  was 
deleted  and  regressions  analysis  done.  The  plots  and  correlation  coefficient 
for  3  minimized  surface  effect  on  the  granulated  powders  is  given  in  Figures 
13-15  while  the  plots  and  correlation  coefficients  to  show  the  miniscule 
surface  effect  on  the  spheres  is  given  in  Figures  16-18.  As  can  be  seen  from 
Figures  13-15  the  pseudo-kinetic  rate  expression  takes  on  the  form  Ca  =  Kg 
exp  (kt)  with  good  correlation  when  the  surface  effect  minimized.  Several 
equations  were  tried  in  order  to  find  the  apparent  kinetics  for  the  surface 
area  effect  for  the  granulated  powders  and  thus  generate  a  predictive  curve. 
The  equation  having  the  highest  consistent  correlation  was  -\j~Ca  =  Kg  t 
The  plots  and  correlation  coefficients  generated  using  the  above  equation  are 

V 

given  in  Figures  19-21.  As  is  seen  from  the  high  correlation  coefficients 
the  above  expression  simulates  the  surface  area  effects  quite  well.  By  now 
combining  the  results  of  the  previous  plots,  it  is  seen  the  spherical  shaped 
aluminas  fluoride  removal  can  be  simulated  by  Ca  =  Kg  exp  (kt)  integrated 


rate  expressions  up  to  constant  removal  rate  region  and  the  granulated  powder 
by  wTC^  =  Ke  t  k  0-2  hours  then  Ca  =  Ke  exp  (kt)  from  2  hours  to  constant 

removal  rate  region  start.  This  is  shown  below  in  Figures  22  and  23  for 

( 

spherical  and  granulated  power  aluminas,  respectively. 

i 

TIGURE  22:  PREDICTIVE  FLUORIDE  ION  CONCENTRATION  IN  ALUMINA  COLUMN  VS 
ELAPSED  TIME  MODEL  i 


\ 


FIGURE  23:  PREDICTIVE  FLUORIDE  ION  CONCENTRATION  IN  ALUMINA  COLUMN  VS 
ELAPSED  TIME  MODEL 


CONCLUSIONS 

The  predictive  models  for  finding  fluoride  concentration  during  the 


Jpputial  start-up  (switch-over  tiroes)  allow  the  alumina  column  process  to  be 
more  closely  monitored.  Using  some  of  the  Ruebel  and  Hager,  Inc.  start-up 
data  further  verification  of  these  predictive  models  was  made.  The  plots  of 

i  .  • 

In  concentration  vs  Time  and  the  correlation  coefficients  are  given  in 
Figures  24  and  25.  NOTE:  In  Ca  =  Ke  +  ktj  is  equivalent  to  Ca  =  Ke  exp  (kt). 
As  is  seen,  the  correlation  coefficients  are  good  considering  the  surface 

l 

area  effect  was  not  considered  in  these  analyses.  These  pseudo-kinetic 

t 

rate  test  results  are  not  completely  conclusive  without  several  repetitions 
of  the  tests  made  on  the  aluminas,  A  further  pseudo-kinetic  rate  test 
coupled  with  column  test  should  be  made  to  enable  fluoride  removal  in 
alumina  columns  to  be  predicted  by  the  mathematical  models  found. 


FIGURE  1 :  REGRESSION  (f  (CA)  US  kt)  on  ALCOA  F-l  (VIRGIN) 

SELECT  BEST  FIT 


, EQUATION .Y  :*  A*XtB  HAS  MINIMUM  MAXIMUM  ABSOLUTE  RESIDUAL 


FIGURE  2:  REGRESSION  (f  (CA)  US  fct)  on  ALCOA  F-l  (REGENERATED) 
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EQUATION  Y  *»  A  +  B*L0G<X>  HAS  MAXIMUM  R-SQUARE 

EQUATION  Y  «  A  +  B*LOG<X>  HAS  MINIMUM  MAXIMUM  ABSOLUTE  RESIDUAL 


FIGURE  3:  REGRESSION  (f  (CA)  US  kt)  on  MCB  SPHERES  (0  X  14) 
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EQUATION  y  =*  A  +  B*X  HAS  MINIMUM  MAXIMUM  ABSOLUTE  RESIDUAL 


FIGURE  6;  REGRESSION  (f  (CA)  US  kt)  on  KAISER  A-300  (UNGRADED) 
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EQUATION  Y  =  A  +  B/X  HAS  MINIMUM  MAXIMUM  ABSOLUTE  RESIDUAL 


FIGURE  7:  CA  vs.  KE  exp  (kt)  -  ALL  POINTS  INCLUDED 


FIGURE  8:  CA  vs.  KE  exp  (kt)  -  ALL  POINTS  INCLUDED. 
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FIGURE  9;  CA  ys\  K.E  exp.  (kt)  «  f\^l  POINTS  INCLUDED 
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FIGURE  10;  CA  vs.  KE  exp  (kt)  -  ALL  POINTS  INCLUDED. 
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FIGURE  16:  SURFACE  AREA  EFFECT  ON  KINETIC  RATE  EXPRESSION  CORRELATION  COEFFICIENT 


FIGURE  17:  SURFACE  AREA  EFFECT  ON  KINETIC  RATE  EXPRESSION  CORRELATION  COEFFICIENT 
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